% Ðàñ÷¸ò à2, à4
%gama0=.5;  % óãîë ðàñêðûâà çåðêàëà
   n0_5=round(atan(tan(gama0)/2)/delta);
            f0_5=Frup(n0_5+Angls/2+1);

   n_1=round(gama0/delta);
            f_1=Frup(n_1+Angls/2+1);

a4=(12+4*f_1-16*f0_5)/3;
a2=f_1-1-a4; 
for n=-n_1:n_1
  ro(n+n_1+1)=R0/n_1*n;
end
clg;
plot(ro(1:2*n_1+1),Frup(Angls/2-n_1+1:Angls/2+n_1+1),'m-');title('ýïþðà');
axis([-R0 R0 0 1]);
%hold on
%plot(-ro(1:n_1),Frup((Angls/2-n_1):(Angls/2-1)),'m-');title('ýïþðà');
%hold off
pause
close
E=f_1/Frup(Angls/2+1)
pause

